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ABSTRACT 


During the Voyager 2 flyby of Neptune, a series of 16 whistler- 
like events were detected by the plasma wave instrument near closest 
approach. These events were observed at radial distances from 1.30 to 
1.99 Rfj and magnetic latitudes from -7° to 33°. The frequencies ranged 
from 6.1 to 12.0 kHz, and the dispersions fit the Eckersley law for 
lightning-generated whistlers. Lightning in the atmosphere of Neptune 
is the only known source of such signals. The frequency range of the 
whistlers (up to 12 kHz) indicates that the local electron densities 
are substantially higher (N g > 30 to 100 cm“3) than indicated by the in 
situ plasma measurements. The dispersion of the whistlers is very 
large, typically 26,000 sec Hz 1/2 . Based on existing plasma density 
models and measurements, the dispersions are too large to be accounted 
for by a single direct path from the lightning source to the 
spacecraft. Therefore, multiple bounces from one hemisphere to the 
other are required. The most likely propagation path probably involves 
a lightning source on the dayside of the planet, with repeated bounces 
through the dense dayside ionosphere at low L-values. 
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I. INTRODUCTION 


The Voyager 2 spacecraft, which flew by Neptune on August 25, 

1989, included a plasma wave instrument that provided the first 
measurements of low frequency (10 Hz to 56 kHz) plasma waves in the 
magnetosphere of Neptune [Gurnett et al., 1989]. During the flyby of 
Neptune, the plasma wave instrument detected a series of highly 
dispersed signals that are believed to be whistlers produced by 
lightning in the atmosphere of Neptune. The purpose of this paper is 
to describe the observations and to discuss the initial interpretation 
of these signals. 

Whistlers, as is well known, are electromagnetic waves excited by 
lightning that propagate through the magnetosphere at frequencies below 
the electron cyclotron frequency and electron plasma frequency. The 
mode of propagation of these signals is known as the whistler mode 
[Stix, 1962]. Since the index of refraction of the magnetospheric 
plasma depends on frequency, the impulsive signal from a lightning 
flash is converted into a whistling tone as the wave propagates through 
the plasma, hence the term "whistler." The wave energy tends to be 
guided along the magnetic field and in some cases can bounce back and 
forth along the magnetic field lines from one hemisphere to the other. 
Terrestrial whistlers were first reported by Barkhausen [1919]. The 
dispersive properties of these signals were later studied in consider- 
able detail by Eckersley [1935], and the first comprehensive theory of 


whistler propagation was developed by Storey [1953]. For a review of 
the observations and theory of terrestrial whistlers, see Helliwell 
[1965]. 

Prior to the Neptune encounter three of the outer planets, 
Jupiter, Saturn and Uranus, were known to have atmospheric lightning. 
The discovery of lightning at these planets came from a combination of 
imaging [Smith et al . , 1979], radio measurements [Kaiser et al., 1983; 
Zarka and Pedersen, 1986], and whistler observations [Scarf et al., 
1979; Gurnett et al., 1979; Menietti and Gurnett, 1980; and Kurth et 
al., 1985]. Since lightning appears to be a common feature of the 
outer planets, it is not surprising that lightning has now been 
discovered at Neptune. In fact, Borucki [1989] has made specific 
predictions regarding the occurrence of lightning in the atmosphere of 
Neptune. According to his calculations, the lightning flash rate at 
Neptune is estimated to be 2.0 x 10~ 4 flashes km -2 yr _1 , which is a 
factor of 19 smaller than at Jupiter. 
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II. OBSERVATIONS 


During the Voyager 2 flyby of Neptune, a total of 16 whistler-like 
events were detected by the plasma wave instrument. These events were 
all identified using data from the wideband waveform receiver. The 
wideband receiver provides measurements of the electric field waveform 
over a frequency range from 50 Hz to 12 kHz at a sample rate of 28,800 
samples per second. For a description of the plasma wave instrument, 
see Scarf and Gurnett [1977]. Since the data rate generated by the 
wideband receiver is very high, 115.2 kbits/s, these measurements can 
only be obtained during a limited number of periods called wideband 
frames. A wideband frame normally consists of 48 seconds of data. 
However, because of time sharing with the imaging system, in some cases 
the duration of the wideband signal at Neptune was as short as 5 
seconds . 

The wideband frames obtained in the region around closest 
approach, which is where the whistler-like signals were observed, are 
summarized in Figure 1. This illustration shows the spacecraft 
trajectory as a function of radial distance, R, and magnetic dipole 
latitude, X m . The radial distance is given in Neptune radii, where 
1 R N = 24,762 km. The magnetic latitude was computed using the offset 
tilted dipole (OTD) magnetic field model described by Ness et al. 
[1989], using the most recent parameters (known as OTD 2) [N. F. Ness, 
personal communication, 1990]. The circles along the trajectory show 


6 


the locations of all the wideband frames obtained near closest 
approach. All of these are 48-second frames, except for the 5 frames 
on the inbound leg near 60-70° magnetic latitude, which are 5-second 
frames. The open circles indicate frames in which no whistlers were 
observed, and the solid black circles indicate frames in which 
whistlers were observed. Each whistler has been given an 
identification number, from 1 to 16. These numbers are listed in 
Figure 1 next to the wideband frame in which the whistler was observed. 
As can be seen, all of the whistlers were observed near the planet, R 
< 2 R n> and at low magnetic latitudes, -7° < X m < 33°. The whistler 
occurrence rate, computed by dividing the total number of whistlers by 
the total amount of wideband time, is approximately 1.6 whistlers per 
minute. 

The whistler-like events listed in Figure 1 were all identified by 
visually inspecting frequency-time spectrograms of the wideband data. 
These spectrograms were generated by Fourier-transforming successive 
60-msec segments of the electric field waveform, and displaying the 
Fourier amplitudes in the form of a spectrogram. Any feature, no 
matter how weak, that consisted of a narrowband tone decreasing in 
frequency with increasing time, was identified as a whistler. Roughly 
half of the events are very weak and require special viewing equipment 
(back-lighted slides) to be seen. The remaining events are clear and 
easily identifiable on normal black and white or color prints. 

Two spectrograms with easily identifiable whistlers are shown in 
Figures 2 and 3. Figure 2 shows whistler number 4, which occurred at a 
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radial distance of R = 1.44 R N and a magnetic latitude of X m = 21°. 

This is the best example of a whistler-like signal detected during the 
entire flyby. Even then the signal is very weak. The maximum electric 
field strength of the whistler in Figure 2 is estimated to be only 12 
pV/m. The narrowband character of the signal, decreasing monotonically 
in frequency with increasing time, is clearly evident. The frequency- 
time slope of the signal, df/dt, also decreases with increasing time, 
as is characteristic of low frequency whistlers in the Earth's 
magnetosphere (see Helliwell [1965]), and in Jupiter's magnetosphere 
(see Gurnett et al . [1979]). Figure 3 shows whistlers number 5 and 7. 
These signals occur in the same general frequency range as whistler 
number 4, but over a somewhat narrower frequency range. Whistlers 
number 6 and 8 also occur on this same spectrogram, but are too weak to 
be seen in this display. 

Eckersley [1935] showed that the arrival time, t, of a whistler at 
a frequency, f, is given by a simple relation called the Eckersley law, 

t = D/VF+ t Q (1) 

where D is a constant called the dispersion, and t Q is the time of the 
lightning flash. To provide further evidence that the dispersive tones 
detected by Voyager are in fact whistlers, we have measured the arrival 
time of each whistler as a function of frequency and fit the Eckersley 
law to these arrival times. Of the original 16 whistlers, only 8 could 
be measured with sufficient accuracy to provide reliable least-mean- 
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square fits. In general, the fits to the Eckersley law are quite good, 
although in some cases the frequency range over which measurements 
could be made is quite limited. The best fit D and t Q values are 
listed in Table 1. To provide an overview of the measured dispersion 
characteristics, Figure 4 shows the arrival time (dots) for all of the 
whistlers for which adequate fits could be obtained, plotted on a 
common time base, t - t Q . Also shown are some representative 
dispersion curves (solid lines) computed from the Eckersley law. The 
results are quite striking. All except whistler number 1 are seen to 
be bunched together along a narrow range of dispersion curves, with D 
values ranging from about 24,300 to 29,200 sec Hz 1/2 . Furthermore, 
even though the frequency range of some of the individual events is too 
narrow to provide a good test of the 1/vTlaw, the overall fit is very 
convincing. In particular, the slope, df/dt, of the individual events 
varies with frequency in almost precisely the manner predicted by the 
Eckersley law. This good agreement provides convincing evidence that 
the signals are produced by dispersive whistler-mode propagation from 
an impulsive source. As far as we know, lightning in the atmosphere of 
Neptune is the only plausible source of such signals. 
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III. DISPERSION ANALYSIS 


Since the propagation path and dispersion of a whistler depends on 
the magnetic field topology and plasma density distribution, whistlers 
can provide useful information on the structure of Neptune's 
magnetosphere. At low frequencies it can be shown [Storey, 1953] that 
the ray path of whistler mode waves is confined to within 19° of the 
magnetic field direction. Therefore, to a first approximation the ray 
paths follow a magnetic field line. It can also be shown [Stix, 1962] 
that the whistler mode cannot propagate at frequencies above the 
electron cyclotron frequency, f c = 28 B (Hz), where B is the magnetic 
field strength in nT, or above the electron plasma frequency, fp = 
9000\/N (Hz), where N is the electron density in cm“3. Since whistlers 
were observed at frequencies extending up to about 12 kHz, these 
conditions imply that the absolute minimum magnetic field strength and 
electron density along the ray path are B m ^ n = 430 nT and N m ^ n = 

1.8 cm - 3. The minimum magnetic field condition indicates that the ray 
paths must be confined to the near vicinity of the planet. Although 
the exact magnetic field configuration near the planet is not yet 
known, we can use the magnetic moment of 0.133 Gauss R N ^, given by Ness 
et al. [1989], to estimate that the ray paths cannot extend to radial 
distances more than about 3 R^- Thus, it is clear that the region 
where the whistlers were observed is characterized by low L-values 
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(most likely L < 3). This conclusion is consistent with the radial 
distances and magnetic latitudes listed in Table 1. 

The dispersion of a whistler is controlled by several factors. 

For the relatively simple Eckersley law to apply, several very 
restrictive assumptions must be satisfied (see Helliwell [1965]). 

These include propagation at relatively small angles to the magnetic 
field, wave frequencies substantially less than the electron cyclotron 
frequency (f << f c ) and electron plasma frequency (f << f ) f and 
relatively high plasma densities (f p 2 >> ff c ). When these conditions 
are satisfied the dispersion constant D in the Eckersley law is given 
by the following integral evaluated along the ray path. 




(2) 


Although some uncertainty exists in the ray path and cyclotron 
frequency because of uncertainties in the magnetic field geometry, by 
far the largest uncertainty is in the plasma density, which must be 
known to compute f p . Therefore, before attempting to evaluate the 
dispersion, we must discuss the plasma density. 

Relatively little is known about the plasma density in the 
vicinity of Neptune. For reference a detailed discussion of the plasma 
density in the magnetosphere and ionosphere of Neptune is given in the 
Appendix. Basically three types of plasma density measurements are 
available from Voyager 2, (1) in situ plasma measurements, (2) various 
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limits imposed by plasma wave cutoffs and resonances, and (3) radio 
occultation measurements. In the region near closest approach the in 
situ plasma measurements from the Voyager plasma instrument [Belcher et 
al., 1989] give very low plasma densities, typically 10"^ to 10"^ cm - ^, 
with a maximum of 1.H cm“3. These densities are in direct conflict 
with the whistler observations, since they give electron plasma 
frequencies below the whistler frequencies, which is impossible. In 
fact, since the whistlers show little or no deviation from the 
Eckersley law, the electron plasma frequency must be well above the 
upper frequency limit of the whistlers over the entire region where the 
whistlers are observed, otherwise the signals would have a frequency- 
time form known as a "nose whistler" [Helliwell, 1965]. A conservative 
estimate, based on the condition fp^ >> ff c indicates that the electron 
plasma frequency must be at least 50 to 100 kHz, which corresponds to 
electron density of at least 30 to 100 cm“3. The most likely 
explanation of the disagreement with the plasma measurements is that 
the plasma instrument is simply not able to detect the cold, dense 
plasma that exists near closest approach (see the discussion in the 
Appendix) . 

Unfortunately, the plasma wave cutoffs and resonances observed 
near closest approach are difficult to analyze and do not provide a 
clear determination of the electron density. However, they do not 
preclude electron densities of 100 cm“3, 0 r more. For further details 
see the discussion in the Appendix. Radio occultation measurements 
[Tyler et al., 1989], which are characteristic of the ionosphere around 
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the limb of the planet, clearly show that electron densities as high 

2x10^ cm~3 exist in the ionosphere at an altitude of 1400 km, 

decreasing to several hundred cm'3 a t altitudes of 5000 km (R = 1.2 

R N ). The radio occultation measurements also show that the ionosphere 

is very cold. The temperature estimated from the scale height is only 

950° ± 160°K. On the dayside of Neptune the ionospheric electron 

densities are likely to be even higher. Models of Neptune's ionosphere 

[Shinagawa and Waite, 1989] give maximum electron densities of 10^ cm" 

3 

, and comparisons with Uranus indicate that the density could be as 
high as 6x10 5 cm'3 (f or further details see the discussion in the 
Appendix) . 

Armed with these overall impressions of the electron density in 
the vicinity of Neptune, we can now return to the evaluation of the 
whistler dispersion, given by Equation (2). To obtain a rough first- 
order estimate we compute the path length, f, that would be required to 
obtain a given dispersion, assuming that the electron cyclotron 

frequency and plasma frequency are constant. The required path length 
is given by 


vr 

( = 2c — — D . (3) 


To illustrate the possible extremes, we consider two propagation paths, 
one through the magnetosphere, where the electron density is relatively 
low, and the other through the ionosphere, where the electron density 
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is relatively high. The assumed parameters are listed in Table 2, 
along with the corresponding path lengths in Neptune radii ( 1 Rjj = 
24,762 km). For the dispersion we assume D = 26,000 sec Hz 1 ^, which 
is typical of the whistlers listed in Table 1. For the magnetospheric 
propagation path we use f c = 28 kHz (B = 1,000 nT), which would 
correspond to a propagation path extending out to about 2.4 Rjj from the 
center of the dipole. This propagation path gives roughly the lowest 
cyclotron frequency that would be consistent (f << f c ) with the 
Eckersley law. For the electron density two cases are given. The low 
density case (N = 30 cm"3) corresponds to roughly the lowest electron 
plasma frequency that would be consistent (fp^ >> ff c ) with the 
Eckersley law. The high density case represents a relatively high 
density (10^ cm“3) that is comparable to the Io plasma torus at 
Jupiter. For the ionospheric propagation path we have used the plasma 
frequency, 7 MHz, reported by Zarka and Pedersen [1986] for the dayside 
ionosphere of Uranus. This value is probably a little high for Neptune 
(because of the reduced solar UV flux), but not implausible. For the 
magnetic field strength two cases are given, one (B = 10,000 nT) 
corresponding to the lowest field estimated by Ness et al. [1989] at 
the planet's surface, due to the offset of the dipole, and the other 
(B = 100,000 nT) corresponding to the highest field at the planet's 
surface. 

As can be seen from Table 2, the required path lengths are very 
large. They range from 2,159 Rjg for the low density magnetospheric 
path, to 48 R N for the weak field ionospheric path. These long path 


lengths either imply much larger plasma densities than are indicated by 
any of the presently available models or measurements, or a large 
number of bounces from one hemisphere to the other, or some combination 
of these effects. From strictly a path length point of view, the 
ionosphere provides the most favorable propagation path for explaining 
the large dispersions with the minimum number of bounces. Although the 
path lengths are very long, they are not ruled out, as far as we know, 
by any fundamental consideration. Long whistler mode propagation 
paths, in extreme cases involving as many as 200 bounces from one 
hemisphere to the other, have been observed in the Earth's 
magnetosphere (see Figure 4-8 of Helliwell [1965]). The main 
requirement for a long propagation path is that the plasma temperature 
must be very low so that Landau damping is negligible. 


15 


IV. PROPAGATION PATH 


Since a detailed model of Neptune's magnetic field is not yet 
available, and very little is known about the electron density 
distribution around Neptune, no attempts have been made to carry out 
whistler ray tracing calculations. However, certain general 
characteristics of the ray paths can be addressed. From the previous 
dispersion analysis it appears almost certain that many bounces are 
required from one hemisphere to the other. The reflection at the end 
of each bounce can be due to either the discontinuity in the index of 
refraction at the base of the ionosphere, or the refraction of the wave 
normal angle through 90° at frequencies below the lower hybrid 
resonance (LHR) frequency. This latter type of reflection is 
responsible for what are called "magnetospherically reflected" 
whistlers [Smith and Angerami, 1968], and could only occur in regions 
where the surface magnetic field strength exceeds about 18,500 nT 
(which is possible), so that f LHR > 12 kHz. In the Earth's 
magnetosphere multiple bounce (hop) whistlers are frequently observed 
in which the ray path retraces the same magnetic field line, as in 
Figure 5. This type of field-aligned propagation is made possible by 
field-aligned density irregularities that act to duct the signal along 
the magnetic field line [Helliwell, 1965]. Such ducted field-aligned 
propagation produces a train of whistlers, each with a progressively 
larger dispersion, such as illustrated at the bottom of Figure 5. 
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Since there is no evidence of a train of whistlers at Neptune it is 
unlikely that the signals are bouncing from one hemisphere to the other 
along a single field-aligned duct. Some mechanism must exist for the 
ray path to shift across the magnetic field lines on successive bounces 
so that only a single whistler is observed at the spacecraft. 

Non-ducted whistler propagation is frequently observed in the 
Earth's magnetosphere [Smith and Angerami, 1968]. For non-ducted 
propagation the whistler-mode energy can propagate at a substantial 
angle to the magnetic field direction (up to 19°). The resulting ray 
path can then deviate significantly from the magnetic field lines. If 
the wave vector is in the meridian plane, the ray path tends to shift 
to higher and higher L-values on successive bounces, as illustrated in 
Figure 6. This pattern is characteristic of magnetospherically 
reflected whistlers. For such ray paths the dispersion should decrease 
as the spacecraft moves to lower and lower L-values. Since no such 
systematic trend was observed, it is unlikely that a simple meridional 
ray path geometry can explain the whistlers observed by Voyager at 
Neptune. More likely, the ray paths migrate in longitude, as 
illustrated in Figure 7. This type of longitudinal migration would 
allow the ray paths to make repeated passes through the dense dayside 
ionosphere before being detected on the nightside of Neptune, where the 
density is probably much lower. Even though non-ducted whistlers do 
not retrace the same ray path on successive bounces, they can still 
occur in trains [Edgar, 1976]. Non-ducted whistler trains occur 
because the ray path from the lightning source to the spacecraft has 
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multiple solutions. Whether non-ducted propagation can explain the 
single isolated whistlers observed at Neptune will have to be explored 
by computer ray tracing calculations. 

Finally, we should point out that the magnetic field geometry may 
also play an important role in producing the large dispersion. Since 
the largest electron densities probably occur in the ionosphere, any 
magnetic field configuration that tends to maintain nearly horizontal 
magnetic field lines through the ionosphere will result in ray paths 
with large dispersions. The offset of the magnetic dipole from the 
center of the planet may help produce this type of configuration. 

These and other ray path considerations will have to be explored with 
detailed ray tracing calculations when a suitable magnetic field model 


becomes available. 
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V. CONCLUSIONS 


In this paper we have presented evidence for the existence of 
whistlers in Neptune's magnetosphere. The observed signals consist of 
a narrowband tone that decreases monotonically with increasing time. 
Sixteen such events were observed, all near closest approach at low 
magnetic latitudes. The frequency-time characteristics of these 
signals fit the Eckersley law, which provides strong evidence that the 
signals are produced by lightning in the atmosphere of Neptune. The 
upper frequency limit of the whistlers, 12 kHz, is in direct conflict 
with the Voyager 2 plasma measurements, which gives electron densities 
too low for the whistler mode to propagate. Most likely the plasma 
instrument is simply not able to detect the cold, dense plasma that 
would be required for whistler mode propagation in this region. The 
dispersion of the whistlers is also very large (D = 26,000 sec Hz 1 ' ), 
which indicates substantial plasma densities somewhere along the 
propagation path, and/or a large number of bounces from one hemisphere 
to the other. The most likely propagation path probably involves a 
lightning source on the dayside of the planet, with repeated bounces 
through the dense dayside ionosphere at low L-values to accumulate the 
required dispersion. Detailed computer ray tracing calculations using 
suitable magnetic field and plasma density models will be required to 
investigate the ray paths involved. The signals are also very weak, 
much weaker than the whistlers observed in Jupiter's magnetosphere 
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(after taking into account the radial distances involved). The low 
intensities, or possibly other geometric considerations may explain why 
no lightning was detected by either the imaging [Smith et al., 1989] or 
radio astronomy [Warwick et al., 1989] instruments. 
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APPENDIX 

Since the plasma density plays a crucial role in the propagation 
of whistlers, it is useful to briefly summarize what is known about the 
plasma density in Neptune's magnetosphere and ionosphere. 

Plasma Measurements . The only direct in situ plasma density 
measurements are from the Voyager 2 plasma instrument [Belcher et al., 
1989]. In general the plasma densities reported from the plasma 
instrument are very low, typically 10~3 to 10“^ cm - ^. The maximum 
plasma density was either 0.6 or 1.4 cm“3, depending on whether H + or 
N + is assumed to be the dominant ion. These maximum values were 
obtained at 0420 SCET, in the same region where the whistlers were 
observed. All of the plasma density measurements obtained in the 
region where whistlers were observed are in direct conflict with the 
whistler observations. The reason is that the electron plasma 
frequency computed from the plasma density measurements is below the 
maximum frequency of the whistlers, which cannot occur for the whistler 
mode of propagation. The most likely explanation of this discrepancy 
is that the plasma instrument was not measuring the total plasma in the 
region near closest approach. During the pass over the polar region, 
the plasma sensor was oriented in such a way as to optimize the 
detection of precipitated auroral particles. In this orientation, the 
sensor is relatively insensitive to plasma arriving from the ram 
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direction, particularly if the plasma is very cold [personal 
communication, J. Belcher, 1990]. 

Plasma Densiti es From Plasma Wave Resonances and Cutoffs . It is 
well known that various cutoffs and resonances in the plasma wave 
spectrum can be used to determine the plasma density. In fact, 
electrostatic emissions at the upper hybrid resonance (UHR) have been 
found to be in good agreement with the in situ plasma measurements near 
the inbound and outbound equator crossings at 10 Rjj and 11 R^j, 
respectively [Gurnett et al., 1983; Barbosa et al., 1990]. 
Unfortunately, in the region around closest approach the plasma wave 
spectrum provides very little definitive information on the plasma 
density. A spectrogram of the plasma wave electric field intensities 
in this region is shown in the bottom panel of Figure 8. This 
spectrogram is essentially an enlarged version of the high frequency 
part of Figure 5 from Gurnett et al. [1989]. No obvious resonances, 
such as electron plasma oscillations or upper hybrid emissions, are 
present that can provide a direct indication of the electron density. 

At frequencies from 10 to 50 kHz the spectrum is dominated by radio 
emissions. For a discussion of these radio emissions, see Kurth et al . 
[1990]. On the inbound pass the radio emission at about 50 kHz has a 
sharply defined cutoff at about 0325 SCET, slightly above the electron 
cyclotron frequency. This cutoff has been identified by Kurth et al . 
as the f R=Q cutoff associated with the free space R-X mode. The fact 
that the cutoff is very close to the electron cyclotron frequency 
indicates that the electron plasma frequency is well below 50 kHz at 
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this point. A second cutoff can be seen below the electron cyclotron 
frequency from 0300 to 0345 SCET at about 10 kHz. At 0345 SCET the 
cutoff increases rapidly, reaching 56 kHz at about 0352 SCET. This 
cutoff could be caused by the low frequency cutoff of the free space 
L-0 mode, or the f^_Q cutoff of the Z-mode. In either case the cutoff 
only gives an upper limit to the electron plasma frequency, since there 
is no guarantee that waves actually reach the cutoff. Of these, the 
Z-mode gives the upper limit to the plasma frequency. This limit is 
given by f p (Max) = [f L=0 (f L= o + f c ) ]1/2 > where f L=0 is the low 
frequency cutoff of the noise band. The solid line labelled f (Max) 
shows this upper limit. The electron plasma frequency must lie in the 
shaded region below f p (Max). 

In the crucial time period from 0405 to 0425 SCET, where most of 
the whistlers were observed, the high frequency radio emissions are 
almost entirely absent. This could indicate either that the electron 
plasma frequency is quite high, much greater than 50 kHz, which would 
block the radio emission via the f L=0 cutoff, or that the spacecraft 
has dropped below the propagation horizon of a remote radio source. 
Although no radio emissions are present in this region, a weak plasma 
wave emission can be seen extending from 1 kHz to as high as 50 kHz. 
Since this noise occurs in the same frequency range as the whistlers, 
it must be propagating in the whistler mode. Most likely it is similar 
to plasmaspheric hiss in the Earth's magnetosphere [Thorne et al., 

1973]. Since the whistler mode cannot propagate at frequencies above 
the electron plasma frequency, this noise places a lower limit on the 
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plasma frequency. This lower limit, labelled fp(Min), is shown in 
Figure 5. As can be seen, this limit indicates that the electron 
plasma frequency probably extended well above 50 kHz in the region near 
closest approach. At the magnetic equator crossing (0420 SCET) Sawyer 
et al. [1990] have reported a narrow, well-defined burst of waves at 
frequencies from 104 to 212 kHz that they identify as electrostatic 
(n + 1/2)f c waves. If these waves are associated with the upper 
hybrid resonance, as is often the case, then the electron plasma 
frequency is probably 200 kHz, or more, in this region. 

On the outbound pass the high frequency radio emission again 
reappears above the electron cyclotron frequency after about 0500 SCET. 
A continuation of this radio emission appears to extend below the 
electron cyclotron frequency from about 0423 to 0515 SCET. Although 
this noise has intensities comparable to the high frequency radio 
emission, which suggests that it is free space L-0 mode radiation, we 
believe that this noise is Z-mode radiation [Kurth et al . , 1990]. The 
reason is that in the wideband frame at 0424:47 SCET a whistler can be 
seen extending above the low frequency cutoff of this noise. This 
situation cannot occur for the free space L-0 mode, since the electron 
plasma frequency, which is the low frequency cutoff of the L-0 mode, 
would have to be below the frequency of the whistler. Therefore, the 
noise must be Z-mode radiation. This interpretation is further 
supported by the fact that the noise has an upper cutoff that follows 
the electron cyclotron frequency (see Figure 8), very similar to Z-mode 
radiation in the Earth's magnetosphere [Gurnett et al . , 1983]. Based 



on this interpretation, the low frequency cutoff of the noise band can 
then be used to place an upper limit on the electron plasma frequency, 
following the same procedure as on the inbound pass. This limit is 
shown by the solid line labelled fp(Max), on the right-hand side of 
Figure 8. From the limits given in the top panel, we conclude that the 
electron plasma frequency started from somewhat less than 10 kHz on the 
inbound pass, increased to well over 50 kHz (possibly as high as 200 
kHz, or more) near closest approach, and then dropped to less than 10 
kHz on the outbound pass. 

Radio Oceultation Measurements . The third method of measuring the 
plasma density in Neptune's magnetosphere is provided by the Voyager 
radio oceultation experiment [Tyler et al., 1989]. These measurements 
were obtained as the radio signal from the spacecraft to Earth was 
occulted by the limb of the planet. By analyzing the phase shift of 
the received signal, a profile of the ionospheric electron density can 
be obtained as a function of altitude. This profile has a maximum 
electron density of N ma x = 2 x 10^ cm"^ at an altitude of 1*100 km, and 
a scale height of about 1800 km. The scale height indicates that the 
plasma is extremely cold. If H+ is the dominant ion the temperature is 
950° ± 160° K. Electron densities on the order of several hundred 
cm~3 were detectable at altitudes up to 5000 km (R = 1.2 Rjj). These 
electron densities are considerably larger than those detected by the 
plasma instrument at a comparable radial distance. 

Since the ray path of the radio oceultation signal passed through 
the ionosphere near the terminator, where the solar illumination angle 
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is very low, these densities are not necessarily characteristic of the 
highest electron densities that could exist in the ionosphere of 
Neptune. The highest densities would be expected on the dayside of 
Neptune, near the subsolar point, where the incident solar illumination 
is the strongest. The best estimate of the dayside ionospheric 
electron densities is from Uranus, where the low frequency cutoff of 
radio signals from lightning [Zarka and Pedersen, 1986] can be used to 
measure the ionospheric electron density. The dayside electron 
densities inferred at Uranus using this technique give N max = 6 x 10^ 

O 

cm J , which corresponds to an electron plasma frequency of f p = 7 MHz. 
Since the atmospheres of Neptune and Uranus are expected to be quite 
similar, and the difference in the solar ultraviolet flux is not large, 
one would expect the dayside electron density to be somewhat comparable 
at Neptune. In fact, recent models of the ionosphere of Neptune 
[Shinagawa and Waite, 1989] give a maximum electron density of N = 

Ilia A 

c o 

10-' cm J , which corresponds to an electron plasma frequency of f = 2.8 


MHz. 
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TABLE 1 

NEPTUNE WHISTLERS 


Frequency 


No. 

Time 

E1R n 1 

1.30 

Lat . 

Mag. 
Lat . 

D(sHz 1/2 ) 

to 

Range 

(kHz) 

1 

0405:35 

61° 

33° 

45,462 

0356:59 

7.4- 7.8 

2 

0406:05 

1.31 

59° 

32° 

-- 

— 

9.3-10.2 

3 

0408:22 

1.37 

54° 

27° 

— 

— 

11.0-11.5 

4 

0410:35 

1 .44 

49° 

21° 

24,748 

0406:43 

8.9-11.6 

5 

0412:50 

1.51 

44° 

16° 

26,361 

0408:30 

9.7-10.4 

6 

0412:57 

1.51 

44° 

16° 

— 

— 

10.6-1 1 .2 

7 

0413:07 

1.51 

44° 

16° 

29,207 

0408:40 

10.8-11.9 

8 

0413:27 

1.52 

43° 

14° 

— 

— 

11.3-12.0 

9 

0415:49 

1.62 

39° 

10° 

— 

— 

11.4-12.0 

10 

0415:54 

1.62 

39° 

10° 

— 

— 

10.8-12.0 

11 

0417:42 

1.68 

36° 

6° 

26,103 

0413:43 

10.3-12.0 

12 

0422:26 

1.87 

30° 

-2° 

27,531 

0417:20 

7.4- 8.2 

13 

0422:47 

1.88 

29° 

-3° 

29,209 

0418:22 

10.0-11.8 

14 

0422:51 

1.88 

29° 

-3° 

24,330 

0417:58 

6.1- 6.9 

15 

0424:51 

1.97 

27° 

-6° 

— 

— 

6.3- 6.8 

16 

0425:30 

1.99 

26° 

-7° 

— 

— 

11.1-12.0 
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TABLE 2 

REQUIRED PATH LENGTHS* 


Path 

N 

B 

f P 

f c 

i 

! 

Magnetospheric 

cm -3 

nT 

kHz 

kHz 

H 

(Low Density) 

30 

10 3 

^9 

28 

2,159 

(High Density) 

10 3 

10 3 

270 

28 

392 

Ionospheric 

- 

- 

- 

- 

- 

(Weak Field) 

6 x 10^ 

10* 

7,000 

280 

48 

(Strong Field) 

6x10^ 

105 

7,000 

2,800 

151 


*Using D = 26,000 sec Hz^^ 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


The trajectory of Voyager II in magnetic latitude, A m , 
and radial distance, R, coordinates. The magnetic 
latitude was computed using the OTD 2 magnetic field 
model [Ness et al . , 1989]. The circles indicate wideband 
frames, with and without whistlers. The whistlers are 
numbered from 1 to 16. 

A frequency-time spectrogram showing whistler No. 4. 

This is the the strongest and clearest whistler detected 
during the flyby. Even in this case the signal is very 
weak, only slightly above the receiver noise level. 

A frequency-time spectrogram showing whistlers No. 5 and 

7 . 

A plot of frequency, f, versus time delay, t - t 0 , for 
all of the whistlers for which it was possible to measure 
the time delay. The reference time t Q was calculated 
independently for each whistler using the Eckersley law. 
For ducted field-aligned propagation a lightning flash 
produces a train of whistlers, each with progressively 
larger dispersion. This type of propagation does not 
account for the Neptune whistler observations, since all 
except one of the whistlers have essentially the same 
dispersion. 
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Figure 6 


Figure 7 


Figure 8 


A representative ray path for non-ducted propagation in a 
meridian plane. This type of propagation also does not 
appear to account for the Voyager observations, since the 
dispersion would be expected to decrease with decreasing 
L-value . 

Non-ducted propagation can also allow the ray path to 
migrate in longitude. This type of ray path would allow 
whistlers from a lightning source on the dayside to reach 
the nightside of Neptune, which is where Voyager was 
located when the whistlers were observed. A dayside ray 
path geometry is favored because of the much higher 
electron density (hence dispersion) in the dayside 
ionosphere . 

An analysis of the maximum and minimum electron plasma 
frequency, fp, from observations of radio emissions and 
whistler-mode noise observed near closest approach. The 
electron plasma frequency must lie in the shaded region 
of the top panel. The location of the (n + 1 /2 ) f c 
emissions is from Sawyer et al. [1990]. 








WHISTLER NO. 4 



START TIME: 0410:23, AUGUST 25, DAY 237, 1989 SCET 
R = 1.4 Rm X = 49.6° Xm = 21.9° 
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MERIDIONAL PROPAGATION 



VOYAGER 


Figure 6 
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LONGITUDINAL PROPAGATION 
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8-G90-375 



RADIO EMISSION WHISTLER MODE NOISE RADIO EMISSION 



SCET 

0300 

0330 

0400 

0430 

0500 

R(R n ) 

3.13 

1.82 

1.21 

2.18 

3.52 


52.1 

69.7 

45.6 

-13.5 

-31.2 
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Figure 8 









